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ABSTRACT 

We present new spectroscopic metal abundances for 74 RR Lyrae stars in uj Cen obtained with FLAMES. 
The well-known metallicity spread is visible among the RR Lyrae variables. The metal-intermediate (Mint) RR 
Lyrae stars ([Fe/H] ~ —1.2) are fainter than the bulk of the dominant metal-poor population ([Fe/H] ~ —1.7), 
in good agreement with the corresponding zero-age horizontal branch models with cosmological helium abun- 
dance Y = 0.246. This result conflicts with the hypothesis that the progenitors of the Mint RR Lyrae stars 
correspond to the anomalous blue main-sequence stars, which share a similar metallicity but whose proper- 
ties are currently explained by assuming for them a large helium enhancement. Therefore, in this scenario, 
the coexistence within the cluster of two different populations with similar metallicities ([Fe/H] ~ —1.2) and 
different helium abundances has to be considered. 

Subject headings: techniques: spectroscopic — stars: abundances — stars: Population II — stars: variables: 
RR Lyrae — globular clusters: uj Cen 



1. INTRODUCTION 

The stellar system uj Centauri (NGC 5139) is the most mas- 
sive and luminous globular cluster of the Milky Way, and is 
the only one showing a clear metallicity spread. Extensive 
spectroscopic surveys performed on large samples of giant 
stars (Norris et al. 1996; Suntzeff & Kraft 1996) have found a 
multimodal distribution of heavy elements. Recent photomet- 
ric surveys have revealed the presence of multiple sequences 
in the color-magnitude diagram (CMD) of uj Cen. In partic- 
ular, a discrete structure of the red giant branch (RGB) and 
subgiant branch (SGB) has been evinced (Lee et al. 1999; 
Pancino et al. 2000; Rey et al. 2004; Ferraro et al. 2004; 
Sollima et al. 2005a), indicating a complex formation his- 
tory. In addition to the dominant metal-poor population (MP, 
[Fe/H] ~ -1.7), three metal-intermediate (Mint) components 
(spanning a range of metallicity —1.3 < [Fe/H] < —1.0) and 
an extreme metal-rich population ([Fe/H] ~ -0.6) have been 
identified (Norris & Da Costa 1995; Smith et al. 2000; Pan- 
cino et al. 2002; Vanture et al. 2002). In the faint part of the 
CMD, Anderson (2002) and Bedin et al. (2004) discovered an 
additional blue main sequence (bMS, comprising ~ 30% of 
all MS stars) running parallel to the dominant one. According 
to stellar models with canonical chemical abundances, the lo- 
cation of the observed bMS would suggest a lower metallicity. 
However, Piotto et al. (2005), from a spectroscopic analysis 
of 34 stars belonging to both MS components, found that the 
bMS stars are ~ 0.3 dex more metal-rich than the dominant 
cluster population, and therefore might be associated with one 
of the Mint components identified from the RGBs. The expla- 
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nation currently proposed to account for the anomalous posi- 
tion of the bMS in the CMD is a significant helium enhance- 
ment (by AF ~ 0.10-0.15) for this stellar population (Norris 
2004). 

RR Lyrae stars (RRL), being representatives of the oldest 
populations, are important tracers of the formation history of 
uj Cen. The first sign of peculiarities in the RRL variables 
of uj Cen was found by Freeman & Rodgers (1975), who ob- 
served a wide spread in metallicity in a sample of 25 RRL. 
This result was later confirmed by Butler, Dickens, & Epps 
(1978) and Gratton, Tornambe, & Ortolani (1986). More re- 
cently, Rey et al. (2000) derived metallicities for 131 RRL 
stars from the hk index of the Caby photometric system. 

We present here new metallicities for 74 RRL stars of 
oj Cen, and discuss the implications of helium enhancement 
on the location in the CMD and the pulsational properties 
of the Mint group of RRL stars. This can be useful in con- 
straining the extent to which a spread in helium may or not be 
present in the cluster. 

2. DATA ANALYSIS 

Observations of 74 RRL variables (38 RRab and 36 RRc) 
were collected on February and March 2005 at the VLT/UT2 
at ESO (Cerro Paranal, Chile) with the multi-fiber spectro- 
graph FLAMES/MEDUSA. We used the high-resolution grat- 
ing HR13 with a spectral coverage of 285 A (6120-6405 A) 
and a resolving power R ~ 22500. For each star we ob- 
tained three exposures (3000 s + 3600 s + 3600 s) secured in 
good seeing conditions (FWHM < 0.8"), reaching an average 
S/N ~ 40 per pixel. The one-dimensional spectra were ex- 
tracted using the GIRAFFE pipeline. Thirty fibers were ded- 
icated to sky observations in each exposure. An average sky 
spectrum was obtained and subtracted from the object spectra, 
by taking into account the different fiber transmission. The 
spectra were then continuum normalized with IRAF. 

A set of 26 metal lines was selected from the database of 
Kurucz & Bell (1995) in the spectral range covered by our 
spectra. These lines were identified, whenever possible, on 
each spectrum, and were fit with Gaussian functions, whose 
centers provided the observed wavelength. The integral of 
the difference between the continuum and the line profile pro- 
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vided the equivalent width (EW) of each line. The abundance 
analysis was performed using the latest version of Kurucz 
(1979) model atmospheres 7 and the MOOG line analysis code 
(Sneden 1973) to compute LTE abundances from individual 
EWs. The abundance was derived from the model that best 
reproduced the observed EWs, for assumed values of temper- 
ature, gravity and microturbulence velocity. The errors on the 
abundances were calculated as the standard deviation of the 
abundance relative to all the measured lines. 

Temperatures were derived from V—K colors during the 
phase interval between the beginning and the end of each ob- 
servation, using the V light curves provided by Kaluzny et al. 
(2004). The mean K magnitude of each star was derived from 
the infrared catalog by Sollima et al. (2004) and corrected for 
the phase effect using the template K light curves by Jones et 
al. (1996). 8 As an independent check, we derived the temper- 
ature of each star directly from its spectrum by minimizing the 
trend in the deduced abundances with the excitation potential 
for the observed Fe I and Fe II lines. The two measurements 
agree with each other to within ~ 100 K for most of the RRL 
in our sample. 

The gravity parameter \ogg was then calculated from the 
universal gravitation law. 9 The microturbulence velocity was 
initially set as £ = 2 km s~ , following Fokin & Gillet (1997), 
and then adjusted within a range of about 1 km s -1 by mini- 
mizing the trend in the deduced abundances with EWs for Fe I 
and Fe II lines. 

We performed the spectral analysis independently on each 
exposure adopting the extreme values reached by the param- 
eters during the observation, and assumed the corresponding 
abundances as upper and lower limits. Spectra taken during 
maximum light were rejected because of non-LTE effects due 
to shock wave phenomena. 

3. RESULTS 

The metallicities derived for the 74 observed RRL stars are 
listed in Table 1. The metallicity distribution is shown in 
Fig. G] The KMM test (Ashman, Bird, & Zepf 1994) per- 
formed on our metallicity distribution rejects the single Gaus- 
sian model at the 99.9% confidence level, supporting a bi- 
modal distribution that peaks at [Fe/H] = -1.72 and -1.22 
(MP and Mint groups, respectively) with an estimated com- 
mon covariance 0.02. A detailed comparison with the metal- 
licities measured by Rey et al. (2000) indicates a systematic 
offset A[Fe/H](this work-Rey) = -0.06 dex, with a disper- 
sion of ~ 0.3 dex. Since the accuracy of both sets of results 
is ~ 0.3 dex, the observed dispersion is consistent with the 
above uncertainties. 

The location in the HR diagram of the RRL stars in our 
sample 10 is shown in Fig. [5] As can be seen, the domi- 

7 The effect of a helium overabundance of AY = 0. 1 on the model atmo- 
sphere is to decrease the continuum opacity, thus enhancing the line strength 
by ~ 5% (corresponding to ~ 0.03 dex; Bohm-Vitense 1979). Given the 
small impact of such an effect, we adopted model atmospheres with a stan- 
dard helium abundance. 

8 Dereddened (V— K)q colors were computed adopting the reddening and 
extinction coefficients E(B-V) = 0.11 (Lub 2002), Ay =3.1 E(B-V) and A K = 
0.38 E(B-V) (Savage & Mathis 1979). Temperatures were then derived using 
the color-temperature conversions by Montegriffo et al. (1998). 

9 The gravity was calculated using the distance modulus (m-M)o = 13.70 
(Bellazzini et al. 2004), the reddening and extinction coefficients listed above 
and the bolometric corrections by Montegriffo et al. (1998). The mass of the 
RRL stars was calculated from the stellar pulsation theory using the relations 
of Marconi et al. (2003, their eq. 2). 

10 We transformed the Mk magnitudes and (V— K)q colors into luminosities 



nant MP RRL population ([Fe/H] ~ -1.7) has a mean lu- 
minosity (log(L/L )) ~ 1.69, whereas Mint RRL ([Fe/H] ~ 
—1.2) are systematically fainter, with (log(L/L Q )) ~ 1.62. 
This evidence is in agreement with the results by Rey et al. 
(2000). For comparison, we show three zero-age horizontal 
branch (ZAHB) models 11 with [Fe/H] = -1.7 and Y = 0.246, 
[Fe/H] = -1 .2 and Y = 0.246, and [Fe/H] = -1 .2 and Y = 0.35. 
The ZAHB models with Y = 0.246 fit nicely the lower en- 
velopes of the respective MP and Mint RRL stars. Conversely, 
the helium-enhanced (Y = 0.35) Mint model is far too bright 
and does not match the observed Mint RRL stars. 

Fig-EJshows the location of the RRab Lyrae in our sample 
in the period-amplitude diagram. For comparison, the mean 
locus of the M3 (Oosterhoff I type; Cacciari et al. 2005) and 
M92 (Oosterhoff II type) RRab Lyrae variables are overplot- 
ted in this figure. As can be noted, the stars of the Mint and 
MP samples occupy different regions of the diagram: while 
Mint stars do not show significant deviations from the M3 12 
period-amplitude relation, MP stars have longer periods and 
fall nicely on the M92 period-amplitude relation. . In this re- 
spect, the effect of a He enhancement on the Mint RRL would 
be to produce longer periods, as a consequence of the higher 
luminosities involved (Catelan 1996). 

An independent check on the helium abundance can be 
made using the mass-luminosity parameter A (Caputo, Cayrel, 
& Cayrel de Strobel 1983; Sandquist 2000). This parame- 
ter depends on the helium abundance in the sense that, at a 
fixed mass and for a given temperature, increasing the helium 
abundance increases the luminosity, and thereby the period, 
of an RRL star. Although A has a relatively low sensitivity 
to the helium abundance, statistical errors are generally small 
((Ja ~ 0.05), and this method may provide some useful indi- 
cations in regard to (at least) relative helium abundances. We 
derived A values for the 74 RRL stars in our sample using eq. 
8 of Cacciari et al. (2005), and converted these values into 
Y differences with respect to the mean Y of the MP stars as- 
suming 6A/5Y = 1.4 (Sandage 1990; Bencivenni et al. 1991). 
Figure 0] shows the obtained values of AY as a function of 
metallicity. As can be seen, there is no evidence of a system- 
atic trend with metallicity to within the errors (ay = 0.07). 

4. DISCUSSION 

The most important result of this analysis is the confirma- 
tion of the existence of a metal-intermediate ([Fe/H] ~ -1.2) 
RRL population in ui Cen with luminosity and pulsational 
properties that are incompatible with any significant helium 
overabundance. This evidence adds new complexity to the 
scenario proposed to explain the origin of the bMS. In fact, 
while spectroscopic analyses suggest that the bMS (Piotto 
et al. 2005) and the Mint RRL share a similar metallicity 
([Fe/H] ~ -1.2), the results presented in this Letter demon- 
strate that they cannot share the same helium abundance. On 
the other hand, a helium-rich population is not expected to 
produce a sizeable RRL component (Lee et al. 2005), unless a 

and temperatures using a distance modulus (m-M)o = 13.70 (Bellazzini et al. 
2004), as well as color transformations and X'-band bolometric corrections 
by Montegriffo et al. (1998). The colors were corrected for the non-static 
atmosphere using the correction by Bono, Caputo, & Stellingwerf (1995). 

11 The ZAHBs were calculated using the FRANEC code (Straniero, Chi- 
effi, & Limongi 1997). The location in the CMD was interpolated from the 
tracks of stars with masses ranging from 0.52 to 0.80 Mq . 

12 Note that M3 is a good reference cluster for the Mint RRL population 
in spite of its lower metal abundance, because metallicity does not have a 
significant effect on the period-amplitude distribution in globular clusters be- 
longing to the same Oosterhoff group (Cacciari et al. 2005). 
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fairly large age difference with respect to the dominant popu- 
lation of uj Cen ( Af > 4 Gyr, in the sense that the He-enriched 
population should be younger) is assumed, in contrast with the 
results by Hilker et al. (2004) and Sollima et al. (2005b), who 
both found a rather small (< 2 Gyr) age difference between 
these two populations of uj Cen. 

Therefore, following these considerations, two populations 
with similar metallicities but very different helium abun- 
dances seem to coexist within the cluster. This possibility is 
still compatible with the complexity of the population mix ob- 
served at the SGB level (see Ferraro et al. 2004; Sollima et al. 
2005b), and sets severe constraints on the helium enrichment 
mechanism and its timescale. 

As a matter of fact, in the framework of the self-enrichment 
scenario, the huge helium enrichment required to explain the 
anomalous position of the bMS in the CMD has to occur after 
the entire MP and (at least part of) the Mint population (com- 
prising ~ 70% of the whole cluster population) were already 
formed without any significant helium enhancement. Then, 
the Mint He-rich population has to be formed from a gas en- 
riched in helium by AT ~ 0.10-0.15 maintaining its metal 
abundance practically unchanged. Under this assumption, if 
the sources of the helium enrichment are intermediate-mass 
stars with M ~ 3 -6M Q (as suggested by D'Antona & Caloi 
2004), then the entire chemical enrichment of the system has 
to occur over a time-scale of ~ 50 - 300 Myr. Moreover, a 
very efficient mechanism is required in order to homogenize 
and efficiently reuse all the "enriched" material ejected by the 
previous generation of polluting stars (Norris 2004). In this 
respect, further support to this scenario could come from new 
evidence suggesting that a similar helium self-enrichment 
may have occurred also in the metal-homogeneous globular 
cluster NGC 2808 (Lee et al. 2005; D'Antona et al. 2005; 
see also Catelan 2006 for a recent, critical discussion). Also 
in this case, no He-rich RRL have been found in the cluster 
(Corwin et al. 2004; D'Antona et al. 2005) 

Alternatively, scenarios outside the "pure self-enrichment" 
one may also be considered; in particular, one possibility 
could be that the Mint populations formed in a different en- 



vironment, thus not partaking the chemical enrichment pro- 
cess of uj Cen. In this case, a complex interplay of chemical 
and dynamical evolution has to be taken into account, includ- 
ing gas exchange with the Milky Way and/or minor mergers, 
within a framework such as the binary cluster mergers sce- 
nario (Minniti et al. 2004). 

Alternative scenarios which try to reproduce the anomalous 
location of the bMS without invoking a helium overabundance 
are not much simpler and present even more obvious incon- 
sistencies with observations. For example, if we assume that 
Mint RRL and bMS belong to the same population (with the 
same cosmological helium abundance), then the anomalous 
position of the bMS can only be explained by assuming a 
different (longer) distance modulus [A(m-M)o ~ 0.4 mag, 
corresponding to ~ 1.1 kpc behind the main body of the clus- 
ter]. This would be inconsistent with the Mint RRL luminos- 
ity shown in Fig. 2, which suggests that Mint RRL are at the 
same distance as the MP population. Moreover, Piotto et al. 
(2005) found that bMS stars show the same radial velocity 
distribution of the MP group, implying that the bMS and the 
main body of to Cen are at least dynamically bound to each 
other. 

The results presented in this analysis add new questions in 
the interpretation of the complex metal enrichment history of 
uj Cen. However, this additional piece of information is im- 
portant to define a scenario able to explain the inhomogeneous 
chemical enrichment of the stellar populations of this system. 
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FIG. 1 . — Metallicity distribution for the 74 observed RRL stars. 
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FIG. 2. — HR diagram for the target RRL stars of u> Cen. Open circles mark stars with [Fe/H] < -1.35, filled circles mark stars with [Fe/H] > -1.35. Different 
lines indicate the ZAHB luminosity for the indicated metallicity and helium abundance combinations. 
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FIG. 3. — The period- V amplitude relation for the two metallicity groups of RRab Lyrae in ui Cen. The solid line indicates the mean locus of the M3 variables, 
while the dashed line indicates the mean locus of the M92 variables. 
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Fig. 4.— 



Relative He abundances computed via the A parameter for the observed RRL stars as a function of metallicity. Symbols are as in Fig. 2 
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TABLE 1 

Metallicities and properties of the target RRL 



OGLE ID 


[Fe/H] 


""[Fe/H] 


(V) 


<*> 


i°g(r> 


log <L/L Q ) 


RRL type 


A V 


P(d) 


5 


-1.24 


0.11 


14.745 


13.295 


3.804 


1.626 


1 


1.15 


0.515 


11 


-1.61 


0.22 


14.534 


13.158 


3.813 


1.704 


1 


0.65 


0.565 


15 


-1.68 


0.18 


14.397 


12.765 


3.784 


1.781 


1 


0.70 


0.811 


16 


-1.65 


0.46 


14.562 


13.409 


3.839 


1.674 





0.49 


0.330 


20 


-1.52 


0.34 


14.579 


12.953 


3.784 


1.708 


1 


1.12 


0.616 



Note. — 



The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 



